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SirMI'lARY 



Modern aircraft design trends have been toward 
aircraft with long, dense fuselages, smaller aero- 
dynamic surfaces, swept wings and higher wing loadings o 
These aircraft are required to fly at high altitudeso 
Flight data from such aircraft have indicated that the 
linearized dynamic equations do not always predict the 
dynamic behavior accurately. In addition, certain cross- 
coupling effects were developing between the lateral 
and longitudinal motion of the aircraft. 

The effect of the trend toward concentrating the 
mass in the fuselage was investigated by making an 
engineering approximation of the change in the moments 
of inertia caused by this change in mass distribution. 
Four configurations were considered, each one approach- 
ing nearer the limit of a missile or body of revolution. 
Since the lower lift curve slopes and smaller wing areas 
of the swept-wing aircraft sometimes require flight at 
very high angles of attack, the effect of inclination 
of the principal longitudinal axis to the flight path 
was investigated. The effect of sweep was approximated 
by varying the effective dihedral parameter, C. ^ 

In addition the effect of varying , the ^weather- 

cock-stability'* parameter, was studied since long fuse- 
lages affect this parameter adversely. 

The investigation was made using a Reeves Electronic 
Analogue Computer (REAC) by setting up a mathematical 
model of a representative military-type airplane flying 
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at a Mach number of 0,6 at sea level and at 40,000 
feet. The mathematical model was based on equations 
of motion which were extended to include products of 
the perturbation velocity terms. The effect of changing 
the mass distribution, the inclination of the longi- 
tudinal principal axis, and variation of and 

Ci^ was investigated by changing the mathematical 
o 

model as required by the parameter change considered 
and obtaining the REAC response to a fixed or standard 
step-function control input. 

The results indicated that the lateral motion 
coupled into the longitudinal motion and that this 
effect was increased as the fuselage became longer 
and denser. Altitude aggravated this coupling as it 
did any trend toward instability caused by a change in 
the other parameters considered. The rudder caused 
greater coupling than the aileron, A positive inclina- 
tion of the principal longitudinal axis caused an 
increase in dynamic stability and a decrease in coupling 
vdiile a negative inclination had the opposite effect. 

A decrease in produced a marked trend toward in- 

P 

stability and an increase in coupling while a decrease 
in had the opposite effect® 

The investigation required by this thesis was 
carried out in the Computer Center of the University 
of Minnesota® 
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INTRODUCTION 

Fighter aircraft of recent design show a trend 
toward long fuselages where most of the mass is con- 
centrated, smaller aerodynamic surfaces, high rates of 
roll, and increased use of sweepback. Flight data 
from such aircraft indicate that the linear equations 
of motion, which result after several simplifying 
assumptions, do not always predict accurately the 
dynamic behavior» In addition, these same design 
trends have required more emphasis to be placed on the 
consideration of the dynamic stability of an aircraft 
in the design stage# In earlier practice, dynamic 
stability usually took care of itself after the static 
stability requirements were met* 

The effect on coupling and dynamic stability 
caused by the parameters principally affected by these 
design trends was investigated by means of a mathe- 
matical model entered on the Reeves Electronic Analogue 
Computer (REAC)* An attempt was made to use a representa- 
tive or typical aircraft of the fighter type* The 
mathematical model was based on equations of motion 
which were extended to include the products of the 
velocity change usually assumed to be zero* 

The mass distribution change was investigated by 
considering four configurations each more nearly 
approaching a missile from the standpoint of its moment 
of inertia properties* The effect of the inclination 
of the principal longitudinal axes was considered at 
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inclination angles of zero, plus ten, and minus ten 
degrees. The effect of and was studied by- 

changing the mathematical model as required by a change 
in these parameters over a range of values. Altitude 
effects were considered by calculating the dimensional 
derivatives for 40,000 feet, changing the mathematical 
model accordingly, and running the program as it was 
done for the sea level case. 

The effect of a parameter change was observed by 
obtaining successive REAC responses to a standard step 
function control input. Only the parameter investigated 
was changed so that any variation in the time history 
response was attributed to the parameter change. The 
data was presented by superimposing these successive 
time histories. The trend of the effect of the para- 
meter change considered could then be easily observed. 
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DEVELOPMENT OF THE PROBLEM 

The classical approach to airplane dynamics is 
outlined here in a very brief mannero The airplane 
was considered to be a rigid body. Six equations 
were required to define its motion; one force equation 
and one moment equation for each axis. These equations 
included the force and moment components due to aero- 
dynamic and mass reactions together with the force of 
gravity and thrust effects. Second order effects were 
neglected and the equations simplified to a degree 
vAiich permitted hand solution. 

Appendix I develops the equations of motion 
used in this investigation. The classical approach 
has been used in general but certain departures have 
been made. In Appendix II the equations developed in 
Appendix I have been simplified by assuming that the 
airplane initially is in straight and level, steady, 
gliding flight. The thrust effects were allowed to go 
to zero since they are not critically essential in a 
dynamic stability study. 

Derivative terras of order higher than one in the 
series expansion of the expressions for the forces and 
moments were neglected. The effect of these terms is 
known to be secondary, data was available only for the 
first derivative, and inclusion of these terms, even if 
available, would require more components than the REAC’s 
provided. A second assumption generally made is that 
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the product of the perturbation velocities is essenti- 
ally zeroo This yields two independent sets of three 
linear equations* This assumption was not made* The 
equations simplified in Appendix II are presented below; 



The axes used and the sign convention observed is 
presented in Table I and illustrated in Fig* S viiich 
heads Table I* 

Ref* 2 points out that motion of the longitudinal 
symmetric type does not induce motion of the lateral 
assymmetric type. An examination of the equation set 
presented above indicates that there are no terms vdiich 
would couple symmetric motion into the assymmetric 
motion if the assymmetric motion were initially zero. 

Duncan in Ref* 2 also points out that when devia- 
tions are so small that their squares and products can 
be neglected, there will be no coupling of the lateral 
motion into the symmetric motion. When the deviations 
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are not small longitudinal motion vri.ll persist alone, 
but a purely lateral motion v/ill, in general, induce 
longitudinal motion» The motion vfill couple into the 
forward velocity through the term Vjl ; the pitching 
velocity through the terms f rx-I t- j ^ and 

1^»- ; and into the normal velocity, through the 
terra » It can be observed that the coupling terms 

are modified generally by moments and products of inertia,, 
The study of the effect of the trend in mass dis- 
tribution toward a long, dense fuselage was made by 
considering four cases starting with a hypothetical 
fighter airplane of conventional typeo The character- 
istics of this basic airplane are presented in Table IIo 
The fuselage was considered to represent 60 per cent 
of the airplane’s mass. It was taken to have a maximum 
width of five feet and to have an eight-to-one slender- 
ness ratio. The pitching moment of inertia calculated 
using these assumptions and the approximation that, for 
purposes of making an engineering estimate, the fuse- 
lage can be considered a long, slender cylinder resulted 
in a calculated value of the pitching moment of inertia 
vhich was very near the of the airplane. As the 
mass became completely concentrated in the fuselage it 
was assumed that the slenderness ratio would increase 
to twelve to one, the fuselage width vould remain the 
same, and that the pitching and yawing moments of 
inertia could be calculated using the long cylinder 
approximation. In this final form, the aircraft would 
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be a missile vdth X g. equal to and very small. 

Fig, 10 illustrates the engineering approximation 
that was made of the variation in moments of inertia 
as the fuselage became longer and denser. The mass 
distribution changes were made \diile the aerodynamic 
damping terms, control power terms, mass and other 
effects were considered constant. This would, of 
course, be unrealistic but it has the advantage of 
studying the effect of a change in one parameter in 
the "partial derivative" sense. 

The product of inertia term resulted largely from 
the fact that stability axes were used. When the 
% -axis is so oriented in an airplane in the steady 
flight condition that it is parallel to the relative 
wind, the Eulerian axes are referred to as stability 
axes, NACA TR 5^9 indicated that the terms were 

unimportant for conventional aircraft, however, Ref. 7 
showed that the term may have a pronounced effect 
on the lateral stability of some high speed airplanes. 

It has already been noted that the Xyg term is involved 
in the coupling of the lateral and longitudinal motion. 
If the longitudinal axes are inclined at an angle 
to the flight path, the moments of inertia about the 
flight path axis and the axis normal to the flight 
path are: 

I, • ij, 5^'n 
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The angle largely results from the wing incidence 
angle and the angle of attacko 

The trend toward swept wings and long fuselages 
increasingly important » Long fuse- 
lages adversely affect , the "weather-cock" 
parameter. Swept wings have very large values of 
effective dihedral thus affecting, 

effective dihedral parameter. The effect of a change 
in these parameters on coupling and on dynamic stability 
was investigated. The change in the static derivative 
characteristics was selected for study since general 
flying characteristics are influenced primarily by the 
static stability parameters according to Ref. 6. 

The coefficients such as L,j- are dimensional 
stability derivatives. They are direct functions of 
airframe geometry, mass, mass distribution and flight 
conditions; in addition to being direct functions of 
the non-dimensional stability derivatives such as, 

“ These basic non-dimensional stability derivatives 
are, in turn, implicit functions of airframe geometry, 
mass, mass distribution and flight conditions. The 
values of the non-dimensional stability derivatives 
for the basic airplane appear as Table II. The dimen- 
sional stability derivatives are defined in Tables III 
and IV. The values of the dimensional derivatives were 
calculated and appear as Table V. Refs. 1, 2, 3, 7, 

8, and 9 discuss these derivatives, their characteristics, 
and how they are determined. They are assumed to be 
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constant for a certain range of applicability and flight 
conditionso 

The dimensional stability derivatives can be used 
directly as numerical coefficients in the dynamic 
equations of the airframe based on real timeo Thus, 
stability derivatives in this form are used in setting 
up the mathematical model for the analogue computero 
The use of this form of derivative permits the exten- 
sion of the non-dimensional derivatives to other flight 
conditions vdth only secondary error o This was done 
to obtain the mathematical model of the airplane at 
40,000 feet« Certain values of the non-dimensional 
derivatives which were specifically applicable to 
flight at this altitude were used while others were 
extended by the change in the values of ^ and U » 

Phillips in Ref« 4 considered the effect of 
steady rolling on longitudinal and directional stabil- 
ityo The equations >hich he used can be arrived at 
from the equations developed in Appendix I by applying 
the assumptions that he made in order to arrive at an 
equation set which could be handled analyticallyo 
He retained the product terms which couple the longi- 
tudinal and lateral raotionso He found that rapid 
rates of roll caused instability if the directional 
and longitudinal stabilities are different when the 
rolling frequency exceeds the lower of the pitching 
and yawing natural frequencies of the non-rolling 
airplane. The instability existed only while the air- 
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plane was rolling. In airplanes of short span and 
high fuselage density flying at high altitudes, 
this instability may cause dangerous attitude 
changes. Further, he found that vhen an airplane 
rolls about an axis which is not alligned with its 
longitudinal principal axis, forces were introduced 
which tended to swing the fuselage out of line with 
the flight path. These are the inertial forces 
usually neglected vhen the linearized equations are 
used. 

In Ref. 4, Phillips made many assumptions in 
an effort to obtain equations vhi ch were capable of 
hand solution yet would predict the effect on coupling 
of steady rolling. In the investigation conducted 
incident to this thesis, the effects of the three motionSj 
roll, pitch, and yaw, on coupling as produced by step- 
function control inputs were investigated. As many 
influencing factors w'ere retained as the REACTS compon- 
ents would permit. 

Solutions of the equations of motion used in this 
thesis wDuld require the simultaneous solution of six 
differential equations, three of second order, made 
non-linear by product terras such as and pit . 

The labor involved would have been prohibitive. The 
equations were fitted to two REAC’s and four servo- 
multipliers, The solutions were recorded simultane- 
ously on two four-channel Sanborn recorders. 
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EQUIPMENT AND PROCEDURE 

The equipment used to conduct the investigation 
consisted of two REAC's, four servo-multiplier units, 
two four channel Sanborn recorders, and a one-shot 
multi-vibrator time switch* The entire computing 
system is illustrated in the photograph presented as 
Fig, 2, The Sanborn recorder is shown in Fig, 3, the 
lateral motion and control function computer in Fig. 
the longitudinal motion computer in Fig, 5, and the 
time switch in Fig, 6, 

Each of the four servo-mechanisms drove three 
multiplying potentiometers. These multipliers provided 
of the product of the inputs. This restriction 
made it necessary to build up the perturbation velocities 
to ten times their actual value in order to obtain a 
simple product such as pn. , The additional amplifi- 
cation required by the servo-multipliers made computer 
balance very delicate and this was a source of difficulty. 
The Sanborn recorder made time history records of 
four quantities simultaneously. The imprint was made 
by a heated pen on sensitized paper. The point did not 
touch the paper, thus eliminating the drag usually 
associated with a recording pen. 

The time switch was designed to provide an accurately 
timed step-function as a standard input. Basically, the 
switch is a one-shot multi-vibrator vhich fires just 
once when the operating button is pressed. This allows 
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current to flow in the coil of a relay in the plate 
circuit of the "normal ly-off” tube. The relay contacts 
send plus or minus 100 volts through the aileron, rudder 
or elevator function svhLtch, v/hichever is closed. 

The time that the relay remains closed is governed by 
the RC time constant controlled by the variable potentio 
meter. A range of from approximately 0.5 to 2.5 
seconds was available. The electrical schematic is 
presented as Fig. 7® A photograph of the "bread-board" 
model appears as Fig. 6. The contacts of the relay 
carrying plus or minus 100 volts would appear electri- 
cally just ahead of the function switches in the REAC 
solution schematic which appears as Fig. 9» 

The six equations of motion were put in computer 
form in Appendix II. These equations are 11-13 through 
II-IS. A computer schematic for the simultaneous 
solution of these equation was laid-out and appears 
as Fig. 9. From the schematic, patchboards were wired 
for the lateral motion and control computer and the 
longitudinal motion computer. This division of the 
problem was based on the wiring scheme vdiich would use 
the fewest amplifiers. The apparent symmetry was 
incidental. The numerical values of the quantities 
entering the problem were translated into potentiometer 
settings. The values recorded were u , 100s » 
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and lO<^on one recorder; '\f , and Son the 

other recorder. The plus and minus 100 volts was 
placed on one side of the timer relay contacts. The 
other side of the relay contacts was connected to each 
of the control function switches, plus or minus, which- 
ever was required by the schematic. 

The graph representing the estimated mass dis- 
tribution changes as the fuselage becomes longer and 
more dense is presented in Fig. 10. Values of the 
moments of inertia of three configurations were picked 
off this graph. New stability derivatives were cal- 
culated as necessary for these configurations. All of 
the moment derivatives such as L-^- were affected. These 
calculated values were translated into potoitiometer 
settings and appear as Tables VII and VIII. 

The configurations investigated were called 
Case I, II, III and IV. The lateral computer was 
referred to as the "blue” computer and the longitudinal 
computer as the "red" computer. These are occasionally 
referred to in the tables using the abbreviations B and 
R. 

The computers were balanced and the recorders 
calibrated. The patchboards were then installed on 
the computers. A rudder, elevator, or aileron control 
movement was then selected and the recorder tape set in 
motion at the desired speed. Pushing the time switch 
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button sent an accurately timed standard step-function 
through the control switch selected. The solutions of 
the equations of motion in response to a step-function 
input were recorded as time histories of the variables 
on eight recorder channels simultaneously. The next 
case was then entered on the potentiometers and the 
solutions obtained as before. These solutions are 
presented as Figs, 12A, 12B, 13A, 13B, 14A, and I 4 B, 

The dimensional stability derivatives were recal- 
culated using data for 40,000 feet and the required angle 
of attack increase obtained from Table II, Where 
specific non-dimensional derivatives were not available 
for this altitude, those for low altitude were extended 
by using the values of f> and IT at 40,000 feet in the 
formulae of Tables III and IV, These stability deriva- 
tives and other factors of the problem which changed at 
40,000 feet were translated into potentiometer settings 
and run off as before. The settings appear in Tables IX 
and X, The results for Cases I and IV are presented in 
Figs, I 5 A, I 5 B, 16A, 16B, I 7 A and 1?B, It should be 
pointed out that the airplane which was marginal in 
dynamic directional stability at sea-level was unstable 
at 40,000 feet. In order to ”fly” even marginally at 
this altitude the potentiometer was set to 0,0115 
and the potentiometer to zero. This had the effect 
of increasing the ”weather-cock” stability parameter 
about six times and reducing the effective dihedral 
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parameter to zero* The computer procedure was as 
beforeo It should be noted that solutions were 
obtained only for Cases I and IV modified as des- 
cribed above o 

The effect of inclination of the longitudinal 
principal axis to the flight path was investigated 
by calculating values of moments and products of 
inertis for inclination angles of plus ten, zero, 
and minus ten degrees. The equations below obtained 
from Refo 6 were useds 

“(li. -10 /uX,v)con^ 

It was found that only the value varied significantly. 
New values of the coefficients containing Xy^were 
calculated and translated into potentiometer settings. 
These are presented in Table XI. Several sign changes 
were required. These were accomplished by using inverters 
as necessary. The computer was operated as before 
using Case I at low altitude. 

The and programs were accomplished by 

changing these values over the range considered, 
recalculating the affected stability derivatives, and 
translating them into potentiometer settings. This 
was done for Case I at sea level and at 40>000 feet. 

These values are presented in Tables XII and XIII. 

The computer operation was similar to the other programs. 
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RESULTS AND DISCUSSION 

Time histories of the variables u , ur , Q > v .. 

. * • 

<p i 9 f <p , and were recorded as the airplane 
responded to a step function control input for each 
change in a parameter being investigated o To study 
the trend of the effect of the parameter changes, the 
curves were superimposed, each curve being drawn with 
a characteristic lineo These composite solution curves 
are presented as FigSo 12A through 23B» The standard 
step function time plot appears as Figo llo It should 
be noted that the A and B parts of the figures each 
represent one of the four channels which were simultan- 
eously recordedo Where separate lines for each para- 
meter change can not be observed, it should be assumed 
that the curve lies behind the previous curve, or curves, 
and that no significant change was caused by the new 
value of the parameter « 

The mass distribution program is presented in the 

FigSo 12A, 12B, 13A, 13B, 14A, and I 4 B 0 The elevator 

response curves show no significant changes caused by 

the changed moments of inertia® The elevator set-up 

the phugoid motion viiich remained essentially constant® 

The short period motion was virtually dead-beat damped® 

<» 

The amplitude of Q in Case IV is reduced since the air- 
craft was more sluggish to pitch change due to the in- 
creased pitching moment of inertia® A slight dutch- 
roll was induced® This was not anticipated since the 
equations indicated that there would be no coupling 
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from the symmetric motion into the lateral motion if 
there was no initial lateral motion® A possible 
explanation might be that there was some unbalance 
in the servo-multipliers or jitter in the lateral 
motion computer due to small line current variations® 
The same effect might be observed in flight tests 
caused by slightly unsymmetrical control application 
or small trim unbalance® Spiral instability was set- 
up presumably by the same mechanism vAiich caused the 
slight dutch-roll to start® 

The cross-coupling set-up by the aileron and 
rudder motion can be observed in the ur plot® The 
coupling caused by the more assymmetric motion induced 
by the rudder step is several times that of the aileron 
cross-coupling® The aileron* s cross-coupling decreased 
as the values of and Tj. were increased® This was 
unexpected but can possibly be explained by the fact 
that two opposing effects were taking place® The 
coupling increased but the airplane became more slug- 
gish since the pitch and yaw moments of inertia were 
increased® This sluggishness can be observed in the 
progressive decrease in amplitude and frequency of 
and IT o The large decrease in roll moment of inertia 
had very little effect on and ® 

The altitude effects on coupling and dynamic 
stability are presented in Figs® 15A, 15B, 16A, 16B, 
17A, and 17B® Many of the observations made in the 
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sea-level program apply here. The short period 
oscillation damping decreased increasing the period 
of this raotiono The aileron coupling was larger for 
Case IV, that is, as the flis elage became long and dense. 
The same effect was observed in the coupling caused by 
the rudder. This time there was little amplitude 
difference between the aileron and rudder coupling. 

The sluggishness caused by increased and , 
together with decreased aerodynamic damping, made the 
swings progressively larger in amplitude and period. 

The roll angle did not significantly change. 

In the program which considered the effect of 
changing the inclination of the longitudinal principal 
axis to the flight path, the program essentially, 
only the rudder and aileron step input response was 
investigated. As Yj became more positive, dynamic 
stability increased; more negative, dynamic stability 
decreased. At an V| of minus ten degrees, dynamic in- 
stability developed. From the coupling standpoint, the 
same effect was observed. There is less coupling when 
\[ becomes more positive; more when becomes more 
negative. The rudder curves show these effects more 
clearly than the aileron curves. 

The effects caused by varying the ”weather-cock” 
parameter, , is presented in Figs, 20A, 20B, 21A 

r 

and 21B, Only the aileron and rudder forcing functions 
were considered. The curves for two values of 
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at high altitude were superimposed on the sea-level 

curves. At both altitudes similar effects were 

observed. Coupling increased as decreased. The 

dynamic lateral stability decreased as decreased. 

This was to be expected, of course, since decreased 

tail arm, or tail area, lead toward lateral instability. 

As the angular velocity changes increased with this 

instability the coupling increased. 

The curves for the program are presented 

in Figs, 22A, 22B, 23A and 23B, The entire sea-level 

program was presented but only the curve for a 

of 0,20 at 40,000 feet was added. As Cj was increased 

A 

dynamic, lateral instability quickly developed along 
with increasing coupling. The effect was more apparent 
in the motion caused by the rudder forcing function. 

The reduced damping at altitude hastened the develop- 
ment of instability and the coupling increased. The 
curve for a of 0,20 was marginally stable at sea 
level but unstable at 40,000 feet. 

The parameters and had very significant 
effects on dynamic stability. The effective dihedral 
parameter, , was very critical. This would appear 

A 

to make analysis of the dynamics of the airframe very 
important in the design stage of swept-wing and delta- 
wing aircraft. The REAC provides a powerful tool for 
such an analysis. 
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The use of a standard step-function input 
provided a good means for studying the effect of a 
parameter change. However, this method was inadequate 
from the standpoint of bringing out any instabilities 
caused by cross-coupling of the type investigated in 
Ref. 4o 

It should be pointed out that the dynamic be- 
havior observed was that of one airplane , or variation 
of one airplane. Although effort was made to make 
the mathematical model investigated "typical" , it 
would not be reasonable to suggest that all airplanes 
behave in the manner of the airplane investigated. 

On the other hand, the procedures used in the analysis 
were straight-forward and it is felt that the dynamic 
behavior observed might be considered reasonably 
typical. 
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CONCLUSIONS AND RECOMMENDATIONS 

1. Mass distribution changes had no significant 
effect on the motion of the aircraft in the 
longitudinal plane. 

2. The motion of the aircraft in the symmetrical 
plane coupled into the lateral motion in the 
form of a slight dutch roll and also triggered 
the spiral instability mode. However, this is 
believed to have been due to slight assymmetry 
in the computing system. 

3. The assymmetrical motion of the airplane coupled 
into the motion in the symmetrical plane and is 
most significant as a change in or . This would 
be equivalent to a change in the angle of attack. 

4. The coupling effect increased as the fuselage 
became longer and more dense. 

5o The assymmetrical motion caused by the rudder 
had a greater coupling effect than the motion 
caused by the aileron. 

6. At high altitudes both the coupling and the 
trend toward instability increased. 

7. At a positive inclination of the principal 
longitudinal axis to the flight path, there 

was less coupling and increased dynamic stability. 
A negative inclination angle had the opposite 
effect. 
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So A decrease in the ”weather-cock” stability 
parameter, , resulted in a decrease in 

dynamic stability and an increase in couplingo 
9. An increase in the effective dihedral parameter, 

, resulted in a decrease in dynamic 
stability and an increase in coupling. 

10* The standard step-function input which supplied 

the forcing function provided a means for studying 
the effect of a parameter change on the coupling 
and the dynamic stability of the airplane* How- 
ever, it does not appear adaptable to the investi- 
gation of instabilities which may develop from 
the coupling caused by rapid rolling or yawing* 

An investigation might do this by using a pro- 
grammed input which would reverse before the REAC 
limits were reached* 

The dynamic behavior observed in this investigation 
was that of the mathematical model of one airplane or 
variation of one airplane* Effort was made to use a 
mathematical model which could be considered reason- 
ably "typical”* However, it would not be reasonable 
to suggest that all airplanes would respond similarly* 
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SYfIBOLS 

initial airspeed, feet per second 
mass density of air, slugs per cubic foot 
wing span, feet 
wind area, square feet 
weight of airplane, pounds 
mass , slugs 

acceleration of gravity, feet per second per second 
aspect ratio 
Mach number 
wing chord tip, feet 
wing chord root , feet 

moment of inertia about 
axis, slug-feet^ 

Iv moment of inertia about 

*o 

p 

axis, slug-feet 
I- moment of inertia about 
slug-feet^ 

moment of inertia about 
moment of inertia about 
slug-f eet^ 

I, moment of inertia about 
product of inertia with 

axis and axis normal to flight path, slug-feet 



principal longitudinal 
principal tranverse 
principal normal axis, 

2 

flight path axis, slug-feet 
the pitching axis, 

the yaiifing axis, slug-feet^ 
respect to the flight path 



r 

b 

s 

w 

no 

9 

AR 

M 

c. 

Iv 
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p force, pounds 

h moment of momentum, (slug-ft)-ft per second 

C lift coefficient , 

> Cj 3 

Cl rolling moment coefficient — 

° ' Cj sb 

lateral-force coefficient , lateral ferce, 

Cr„ aerodynamic pitching moment coefficient about the c.g. 

Cn drag coefficient , Drag forc^ 

D q s 

angle of bank, radians 
9 angle of pitch, radians 

■\p angle of yaw, radians 

S control deflection, radians 
yS angle of sideslip, radians 

y flight path angle; angle between relative wind 
and the horizontal, radians 
angle between the longitudinal axis of the fuse- 
lage and the relative wind, radians 
^ angle between the thrust line and the longitudinal 
axis of the fuselage, radians 
Ej perpendicular distance from c.g, to the thrust 
line, feet 

y fraction of critical damping 
CO frequency, radians/sec. 

U initial forward velocity in direction of positive 
axis (feet per second) 

Y initial side velocity (feet per second) 

W initial normal velocity (feet per second) 
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u forward velocity change after disturbance (feet 
per second) 

u side velocity change after disturbance (feet per 
second) 

u/ normal velocity change after disturbance (feet 
per second) 

P initial rolling angular velocity (radians per 
seco nd) 

Q initial pitching angular velocity (radians per 
second) 

R initial yawing angular velocity (radians per 
second) 

rolling angular velocity increment, radians per 
second 

» 

C|,0 pitching angular velocity increment, radians 
per second 

yawing angular velocity increment, radians per 
second 

Cr> the change in drag coefficient with varying 
forward velocity 

C the change in drag coefficient with varying angle 

Po< 

of attack, per radian 

rate of change of lateral-force coefficient with 
rolling-angular-velocity factor, per radian 

C, rate of change of lateral-force coefficient with 
yawing-angular- velocity factor , per radian 

Cy lateral force derivative, rate of change of lateral 
force coefficient wdth angle of sideslip, per 
radian 




CL W) 
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C 



L 






c 



L 



oc 



c 



L 




q 







c 



h 



R. 




the change in side force coefficient v/ith variation 
in rudder deflection, per radian 
the change in side force coefficient vdth aileron 
deflection, per radian 

change in lift coefficient vdth variation in forward 
velocity, angle of attack and attitude remaining 
constant 

change in lift coefficient with varying pitching 
velocity, no change of airplane angle of attack 
as a whole , per radian 

the change in lift coefficient v/ith varying angle 
of attack, ”lift curve slope”, per radian 
change in lift coefficient with changes in elevator 
deflection, per radian 

rate of change of rolling moment coefficient with 
yawing angular velocity factor, per radian 
effective-dihedral derivative, rate of change of 
rolling moment coefficient vdth angle of side- 
slip, per radian 

damping in roll derivative, rate of change of rolling 
moment coefficient with rolling-angular- velocity 
factor, per radian 

change in rolling moment coefficient vdth variation 
in rudder deflection, per radian 
aileron effectiveness on aileron power, the change 
in rolling moment coefficient vdth change in 
aileron deflection, per radian 






■fi|aB!OTn^^ 
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C, 



rn^ 



c 






longitudinal static stability derivative, change 
in pitching moment coefficient vd.th varying 
angle of attack, per radian 
the change in pitching moment coefficient \vith 
variation in rate of change of angle of attack, 
per radian 

pitch damping derivative, change in pitching moment 
coefficient v/ith varying pitch velocity, per radian 
Crr,, elevator effectiveness or elevator power, the change 
in pitching moment coefficient \-d.th change in 
elevator deflection, per radian 
directional stability derivative, weather-cock 
stability, rate of change of yawing-moment co- 
efficient with angle of sideslip i radian 

damping in yaw derivative, rate of change of yawing 
moment with yawing angular velocity factor, 
per radian ( YaEt) 

iVo / 

On rate of change of yawing moment coefficient with 

^ \ 
rolling angular velocity factor, per radian ( ^ -p b ) 

* ^ zu„ / 

rudder effectiveness or rudder power, the change 

in yawing moment coefficient with variation in 

rudder deflection, per radian 

Cn change in yawing moment coefficient with change in 

aileron deflection, per radian 

Xy change in fore and aft force with change in forward 

velocity, per second 

Xq change in fore and aft force, change in pitching 
velocity, feet per second-radian 



c 



Hr 
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change in fore and aft force with change in 
normal velocity, per second 

Xj, change in fore and aft force with rate of change 
in normal velocity, per foot 

change in side force with change in yawing velocity, 
per radia n 

change in side force with change in side velocity, 
per second-radian 

Yp change in side force v;ith change in rolling velocityj 
per radian 

Yg change in side force with aileron deflection, 
per second-radian 

change in side force with rudder deflection, per 
second-radian 

change in normal force vdth varying forward velocity, 
per second 

change in normal force with varying pitching velocity, 
feet per second-radian 

Z change in normal force with varying vertical 
velocity, per second 

change in normal force with elevator deflection, 
feet per second^-radian 

[_ change in rolling moment vd.th varying yawing velocity, 
per second-radian 

1_^ change in rolling moment with varying side velocity. 



per second-feet 
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L 



L 






M. 






M 



or 



M. 









V‘ 






M. 












change in rolling moment with varying rolling 
velocity, per second-radian 
change in rolling moment with rudder deflection, 
per second^-radian 

change in rolling moment with aileron deflection, 
per second^-radian 

change in pitching moment with varying forward 

velocity, per second-feet 

change in pitching moment with varying vertical 

velocity, per second-feet 

change in pitching moment with rate of change of 

vertical velocity, per foot 

change in pitching moment with varying pitching 

velocity, per second-radian 

change in pitching velocity with elevator deflection, 
2 

per second -radian 

change in yawing moment with varying side velocity, 
per second-feet 

change in yawing moment with varying yawing 
velocity, per second-radian 
change in yawing moment with varying rolling 
velocity, per second-radian 
change in yawing moment with rudder deflection, 
per second^-radian 

change in yawing moment with aileron deflection, 
per second2-radian 

change in yawing moment with side-slip angle, 



per second'^-radian 




mm 
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( refers to initial steady flight condition 
pertains to the fore and aft axis 
pertains to the transverse axis 
pertains to the normal axis 
C. elevator 
aileron 
( rudder 
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Axfs Designation Longit udincil 


Lotera 1 


VasA/irio 


S y f n bo 1 


X 


Y 


z 


Force ( poraliel b> 

axi's^ synnbof 


X 


Y 
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Sign convention 


^ onward 


To right 4- 


fDo wn 4* 
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linear velocity 
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Sign Convention 
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Charactfristics of the Basic Airpi ane 
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lONGiTUDiNAL Stability Derivatives j 



3 

X 


^(-Co.-Ld) 


X, = 


O 


2u 




2^ 


_e^cc^ 
-4 m 




ei^(c.,cj 

Xa 


Me, 


p USc^’p 

4TT 


Tu 




rpm 


BOpU5c(>^ 

<nn ^Ibrpro 








- PJ^C 

2m ®Se 






2e. 


-pU‘5c^ 

Z nn 




/oU5c r 

2Iy 




Xly '^Sf. 


•II 

X 


0 


Xs. 


-p-^Z5co 

2m 


ZOr 


pSc ^ 
4TJT 


Zs, 


_pU*^Scf~ 
2 m 


Mor 


pSc^r- . 


Ms. 


Z Iy 



Table E? 



Lateral Stability Dbrivatives I 



T.- 


; 


V, 


P USb p 
4m 


Lv 


p-afe., 1 


Lr 


pUSb‘ Q 
4 - lx ^ »- 




Uo L,r 


Nr 


PEt^'Cn. 




p U 5b p- 
2 I. 




P o'-s c. 
Zm 




U.H. 


'-Sp 


2 Ik ^5r 


Tp 




Ns. 


pUT^Cn. 

2 ta Sr 


L., 


p_y^"c, 

4-Tx H 




5 Cy, ' 

2 m c>a. j 


N-p 


pU5b"C_ 

4 - 1 , ^ 


L-s. 


'>!i'’=‘'c,, , 

^ '-S «. 1 






Ns, 


pU*Sbr 1 

ZI> ^‘^aJ 



- 57 



Tob!c m 



DiM£;N 5IONAL DeRlVATiVFb of Th5 B/vsic Airplanl 
M= 0-6 ALTnuoe, Sea Level 



Yu. 1 



-.0151 



Lx-, sec.-fL 



-.054 






see- rad 



L-p) 



I 



Sec - raci 



- 6 55 



X. 



sec 



.0269 



Lg,^ , ssc‘-ro 



rod 



+ 2,84 



Xiir 4-+ 



o 



L 



Sa t sec- rad 



Vr> rad 



+ 1.81 



Y — ^ 

'r> see- 



rad 



M73 



M., sec-^t 



Yluj, sec- f-t 



+-feZ 5 
- 000 ?e2. 
03.:- 






V) rad 



-.302 






I 



^ . sec- roa 



~ 16 66 



Y 



8 ^ » sec- rad 



O 



M 



) sec^- rad 



76 6 



Y 



) scc-rcid 



Z 



u , sec 



*■22 4 



-.‘30 



Z 



-ft 



C), sec -rad 



+“.55 



TT 

I 

sec --ft 

[ 

N se n~ rod 



oo 3d.- 
'*^.o\2S 
~ lOO 



z, 



( 






u/> sec 

4:1 



■ rvd 



L 



P j secured 



-2 70 

+ 61.3 

+- 1735 



Hp. 



see- rad 



I 

^ Aec*- md 



■*■. 0 So 3 
T ^ Z1 
“•5 . *16 



Table V7 



L AT‘^J^Al a^o CoKjr^Oi^ 

^^t^iRur^f^ M/^zz* Dt^Ti^tRuT'iOKj PecOA^ Lc^e^f) 



Pc 




P^^AM7crc^^ 1 






1 


' 


re^rtc^^ere/^ 


^^Rccrtzo 1 


casc r 


C^se iL 


iZL 


CascJZ^ 








~" i 










/ 




0*3 i 


. oos^ 


. 00 2LS 


* Oo f 




2 


r< - 1 *- 

Xf 


' S io 


■ 70^ 


.970 


. ISZ 


I ^ 


— A^q 

■ct<5 


■ 4-3S 


• Z 6 o 


./gs- 


■ 143 


1 

1 

< 

u - 


4 


i- J 


■ g’/f 


■ S>7 


. Elf 


. Elf 


( 

1 _. 


b 




. <^/S* 


. oisr 


. Ofi^ 


. O/i" 






A /3o 


, /9o 


■ ISO 


. !SO 


1 

1 


1 


/ a oo ^o 


.674- 


. <074. 


. 674 


-674 1 




8 


_L_ / 


• 2-97 


■ 291 


■291 


2-77 . 


9 








. (ofZ, 


• 072 . 


10 

u~ 


)(i<r 


. OZ7 


- oz,n 


^ OZl 


. 027 


/ # 


M ur 


*^54- 


.0 24 3 


.OlZis 


. 0 // z 


/ z 


£j, (2.5V 




.OZl 


. 017 


.02-] 


’ /3 






• 0105 


.OII^S 


. 0 1 1'i- 


I 


— A4. 


,oo ^ J 


. 00 0 g 


. OOOfo 


.QOOS 


i t> 


25^ "4;/ 


. 004-7 


A 00 


. 002 ? 


, OOZJ 


/(a 




. ios 


. /OJ 


. /oi 


A /<?5 


11 




* ozo 


A C52o 


. 020 


. OZo 


! 


7^’ 


* ooz s 


* 00 


. 00 4 S' 


-0077 




7^ 


* 0^^ 


. 0*7/ 


. too 


. nz 


ZO 




zs^ 


. 2 5-4, 


. 2S-6 


- Zb4 


1 ~ 

z/ 




. lS(o 


. 2 SC, 


. 2£TC 


‘ZS(o 



- 39 - 

Tahle YIL 

^ itiAj ^err/AJos Fk3j^ LoAj<jir~u0iAjA<~ do P CA y~ e 

Mass C^>zT-f^isu'no'‘j ^p.oc>p.am — deA./! 



fh rvtsj rt erff >v. 


AFFe<.rt o 


CAse r ' 


CAtC: 7L 


CAztm 


CAStJSO 


/ 


Kl-r 


.ot^ 


. 0/0 


■ ooS 


■ ■. — ■■■ : , .. y 

. OOC> 


2 


t^n. 
/ o 


■ ot<i 


i 0/4./ 


, 0 Q n 


. 00<?U 


3 


i i -^> 


. 0^5 


. / 4-Oi' 


./9o 








. 035 


• O / 


. 0097 


. oo-^ 


s 




/ Oo£ 


.COS'] 




.002-^ ; 


(o 


J. ^e- 
4. Tg- 


. ooS 


. 004.^ 


• 00 24 , 


oof 


7 


IX 


. SS^ 


. SB5 


. Si! 


■733 




/ Jvt 
^ lx 


.O!^ 


.01^4 


.0114 


.0/48 




Jxe- 

lx 




>oni 


. ob9S 


,09^15 


to 


!OC> 


.ohn 


.o%i 


,12Z 


.LOf 


U 


"io 


.o/? 


, OZZ-L 


.035 


,o56z 


!Z 


/ ooo 


. In 


.in 


‘ 1/t 


• 1/! 


!3 


fl-<3 


.g”/9 


■ 3t<9 


> f/<7 


\ 

. ^/<f 


/4 




. /S’! 


, /Si 


. t^i. 


> /?2 


/5 




. o3o 


. OSO 


. oio 


• 030 


iCo 


rv 


. Z74. 


. /14 


' 114- 


. 174 


n 


^ V 


.0^4. 


. <3694 


- ofS7 


. /G8 



40 - 



Tc* b Ifc ~V ! IL 

K/ and F'(2.0DUCT5 InJ^RTiA 



Fee. CAf.es J,JL,JJJ 



Case 


r 


TT 


HT 


32: 




5’3 , 000 


4 ^t ^ 000 


Z 0 0 0 


/ 7 d? (i>o 


ly (il ^3 ■ ^^'•) 


2 o, 000 


f 5 4^^0 00 


/ c?oo 


2. C 300 


(slug- ft”) 


/ 2 . 4 -, e 00 


l(ol , 000 


i 

Zn ,000 


2 S’ 0 00 


/ 

Ixt (5I0J 'ft') 


d.^000 


5 , oao 


A 000 

! 


! poo 



lablQ je: 

T T"/ o c? ^ £ 'TT it^ ^ lb 

I, A TC f^, c. 4 J- Co Qo^ t'u TCK^ 

/S tyT “/ C/W P^J 3 u^Af^ Uo, coo rtr-rj 



1 

romt^i-tGf^er^r^ 


i 

C.5. r 1 


C.A^<^ 21L. 


<Z^sa^ III- 


C/iSif. -Zj? 


/ 




— — 1- 

♦ O o tT8 


. OO z8 


A ao f f 


z 






, 9 70 


. 7dz 


3 


,o^ZS 


.05-5^3 


. 0284 . 


^ o3oE 




1 ^ 


, Si^ 


• ?/«? 


■ ?^9 


5 * 


. ^07c> 


, ocjo 


. 0070 


. <50 7fl 


6 


, /:?4L 


> /34^ 


. /W4- 


. /S^ 


7 


.■5 S3 


.^8i 


^ ^6S 


.£83 


s 


. ^3632 


.oCa-SZ 


. 0632. 


.0&3Z 


9 




. t^l4- 


. /^7V 


‘/ 474- 


1.0 


.Ol^Of 


.OZ‘^‘f 


- 02.99 


. Oi99 


// 


. oo-Js, 


. O04-S, 


• 003/ 


,0024- 


12 




. 0099 


- 00 99 


. 00 99 


>3 


.oiZQ 


, OO ~J4- 


. Oo^J 


, oo4og 


i4 


.COOS' 


* Ooo ^ 


. 0005 


.0003 


/^ 


. ooz€ 


.oo/S 


■ 00 14 - 


.OOIZ. 


/^? 


. /03 


,/02 


, 103 


. /OS 


<7 


oo/S 


.00!% 


> 00 /^ 


. 00/3 


•3 


, O OCi 


,000 2. 


. 0003 


. 000 £s 


17 


r OO^ 


p oaj 


. OO^Z 


.0/4-3 


20 




. Z SG 


. Z%(o 


. 2^?6, 


Z! 




. z su, 


. 2^6 


, 1 



5/. - - 



A/o i € ; 



2.5 



- 4.1- 
7^. 31 



c .^yt ^ -S * f Ki j^as 

/^ < ^ 6 frf^'o c Cor^i Pu 

Ma SS ^ # s Tf< f&UT'fO'^ /■ VO C3 /i o r 



Po T» fsj T - 
1 OA^ Brect^ 


ca 5 c r 


c AtbB nr 


Case TTI 


Case j3T 


t 


- Oi IS 


(. oo / 9 J 


^ Oo / 


- 00 S(o 


z 


. oa^l 


. oo 5 


- <30 24- 


■ 00 ?-o 


3 


.OSS' 


. /^s 


!9o 


. Z Z4- 


JL 


. 033 


. O 184. 


.00 91 


. i 


G 


/'” , ' ■ 

/ , / 

/ 


^ / ' / 

' /^ Oa ^4* 


/ / /T 

Ocy^- y 


'/ / ' 

,a.of(o / /' 


4 


zOoS 


- oa dLG 


^ 00 z 4 . 


> Oo ! i 


7 


■ ss^ 


. 6^3 


. 8 Z! 


' 733 






^0/8^ 


. 01 14. 


• 0 / 




.o-]to 


• 07^ 7 


/ <06 <iS 


'OS is 


io 


.o!3>^ 


. °/ 8 ! 


,ozS5 


■ 0 4-S3 


n 


.OBtfS 


. o4-‘^S 


^0(o3S 


/ /057 


'Z. 


. S63. 


’ ^83 


-^83 


. 533 


<5 


■ 8t^ 


. g/f 


.< 5/9 


.811 


id. 


. 0 5S6 


, 0356 


-0386 


■0386 


IS 




/ / op /^ 


/. 00 tS 

/ . . / 


/ . 00 /j 

■' ■ , / / 


/ 4s 


.Oil 


. 057 


.031 


i 

1 

0 


n 


o 


(^0 nfs) 


( .OZ !q) 


0 



in a ic/i tnj io^ic.<xi<^ 
Ghan 



<x St ^ n 




- a-3- 



/ able, ^ I 



a TiO e- 






T x-h f^^oa RAf^ 



n 




^ Tir 

i^tC 


± Ixe- 

~r^ 
8 << 


-X 


L Ixe 
^ '^J 

• B '' 


3 


r% -' Ty 

r< 

_ 7 


Ik -If 
UT 


-lO*" 


. /Oh 


.OZSS 


. O 74 


■ 227 


• oso 


^0 


.g /3 


. ?So 


o 

o 


o 


o 


o 


o 


G 






. 7 oz 


o 

4 - /O 




^O-LS^ 

y y 








X <=>67 


■ S >3 


.SSo 



C/^oz^ hateJ^ s- ^tgn 



/o^/e 221 



C/7 P^OOR^AM 





^TT 0 <) 


>^v (i<) 
4^,000 


. o 


G) 


<G 


t-.o£5‘ 


- <0/3 




/o 


. 023& 


, o<5 73 




.0^54- 




y. 30 


.oJo3 


,0/5^ 



/ ^a^ /e. )LU L 

PRjDCEiR.Ak-\ 

S 3 



“<3 




L-r Ol!^) 




o 


o 


y /<2 


• og-4. 


.OHS 


/2o 


. / o<S 


• OZ3o 


«3o 


^ /c,2. 




/4~o 


I ,2/^ 





XoVo plane 



- 4-4 



Y 






• ^ 

Y 7 



r> 

A 



N Y 



/ 

k, 2 , 



cor? 



X| - ■*■ '^'^■''''P 

Y, * “ XoS»n 4 

^ - 7 

—\ io 



a / -^.c^Q 

V e'^Tr^ 
Yin 6 Y«>ir''V ' 



Yi "Z^ plane 



2 , cos 9 



y 



Z, sn 6 

X-j -Xi Co 5 .&“-i^sin 0 

Yr*-Y 

Zj^ = ?jcos 0 +)(, 5 m 6 




/ 5 . 



Z f 




Xi= C" + CWcosdijSi’rt^o) ^<>^00 ^'I'o) MoS , 

+ Sin-^sm© Sir ^ -V C.Wcas Qa coicib^O^i® ^•<'<^) , 

?i” (rWiiriajK^0i4^ +^>f'|i\n4>) 4CWcQi9a i'r'(^,ystn4' wrS COS<^ 

— cos^* Si -v- O^casGo C05 4) 

Fiq. 1 OtTERMiNATtON OF rrt£ ComPONcNTS P 

v.‘ 

disturbed suuSrjam .\>e 5 



CqRAviTY AcnH<3r AuPNCr 






4 



-k5- 





ThB S4A/^0Aa/ CHA C o 




-U6- 




f,y 4 The care/iAi. At^e> 

COArr^i-oi, Fut^cTtotJ ccif>^f^u reii. 




Ft'S- ^ 'ThC LOi^&trUPtsjAL MOTiOt^l 

O PCJT^O^ 



I 



- 47 - 




Fig. 6 Th-e one- shot Mi/LTf^^/ffAATojz. yir^E St^trcH 



- 48 - 



0 

0 

N 



\n 

vo 




I 






Shot MuLTiv/eRATOfi TJme Switch 



n 

-toi> i 







Z 6 J-- 






2 / 



-7g*^g / 57 A t- 0323 g 



lO 



. 0 S<P 




- 






4 




4 - to 'P 



-A 



Ql 







~ /o*f 



©"^ 



- /O/' 




4 



,© 



0 - 



^/ 39 ^ 






^5 







JO 






©' 

@ 




^ 5 } - 



-/o ^ 



- o SC^TT 










0 * 



^ /O^ 



// 




■ o/r©) 



e-.j 












- <* 




+ /OT> 



© 



/ 




-jX-v-°-* 








^-.U^ ©J" t © 




- /o ^ 




5 u ^>r?er 



In verier 



54 j ^tenif'fo/Tjet’c/* 

se^t/rry 

Computer- No / 
longitudinal G<J *■ control 

Computer No.Z 
lateral ecfuat i ons 

Co rrt rn on points 

©@ 



F /^. 9 REAC Solution 







50 




51 



J o«« 






2**On 






tOOw 


Si, on cJ or~cJi 


pulse ■ O.^Z^ec. 


o 


r— 


L 






o 

i (jec.) 


^ /O 









F I . i . C- . 



U - U I 



I F. •. 



" I 



■ I . . ' 

. I L ^ t.' •' * . ' T U • 



VA ' " ^ C uu.' t . ■ ’ 

■ Cl >< I ^ j ^ > N' • 



t. I C T : A N ' , 

. T F I J U , - L 



K 'I . . T J V . 

i . .. r • * t f I i > •• • 

( I C/^ T : - I . * . ' I ■ 



\ h'-' • • j , : 

T T i ( 



I • t- U ► 'A 
A r< A ‘ 



^ 'O 7 






. 

r . 

( / 



' r f 



C.N'T 



I 



- 52 - 




Itb’ i<fnc Lust 




/o 






6 ^ 












- s 




T 

• 

• 

• 

1 




-/o 


t 

1 








_ _ 














?i5! 

,o2N— 






itb 

O / ^ ^ ^ 






-.02 






-o4 






• 







vA . I. I. I , I 



- 53 - 




O 



/£L 



•>n^ 



"t irrt € 




O ^ / O ‘ ^ 



t C.) 




wo 












11 0 


•Q> 


r. 
















-lo 












1 V_^T' 







r I * 1 ' . ' f'w.C T T ^ V \ - '-< 

' . I , M , I I , A . I y (3 ^ / ' • 



- 54 . 





s “> '-S' 




+ (O 






f- oS 






o 






-05^ 






- . /o 













F I 



V. M \ I ' r P' ^ y 

\ J .. 



9 



1 



V 



55 




• z 






. / 






o 






... ^ 






CD 

-.2 










- 



o if /O '5" 



i C^ec.) 





r I . ; 



1 1 



« I 






3 < 



I 




60 






4-0 /\ 


^ I 




1 

I ^ 1 

Wv \ 


V 7 i' 

\ / ' V 


K\ 






/ i 1 

/ > 1 
/0 1 




7T V 

1/ 




1 A ^ 7 V \ / 

V \ ' A V '/ 


y\\ 


-20 


IU\ 

V\\ 1 




r A 


V \ y \v/> 


\ A^V/ \ 

/\/ 






ra:3Q 


/vy 






f ( sec.) 



fS 



.3 






.2 /\ 






■^1 ^ / /- \ ^ 

//A /A 

.i'/yV- 








"Av^Vn 




i v4" \ 


XA->A\) 




W \\¥^' \ 


vy ^ 





. t$ 






. /o 






0 V 

coJiJ2 






e 

0 






-.o7 













Fi 



< p 

A 






- 57 - 




. z 






f 






'O 

o 

L 






0 






-.2 













O 5 lO >7 





F t ^ ^ TO A ,ri, -• ■ ■< I 0-. , 

AC :: I , I I , 111, Awn I J ' 






Ow A V r y 




i 



! 




- 58 - 





o 


S 


i (sec.J 


fO 


fS 


o 


iOo 




Z act 


^oa 




, ,o:c • r . . 



- 59 - 




lU- 



•on^ itrne base 



to 


















- S' 






~/.o 













o 6" to /iT 

Csec.) 





I 



1 J ^ I F I '3 







r 



- 60 - 





/t 



■i (^ec.j 



to 












0 












-t.o 



















^ 0 
0 ^ 

05 ^ 






•Q> 

O 






-. oS 



















I 



r 



- 61 - 



2 






/ ^ 






:5 

o 




. - 


- 1 


“ . 




-2 













.04 






02. ^ 






03 

O 






~ .02 


m 




-. oJ^ 













O iT to 



t (sec.) 




I A r.'G ^ I r I ' V % , ’ ' . r 



- 62 - 




.2 






■" "V ’ " — ■ ” ■ 

Cf 






o 







-./ 






-.2 













S ,o 'i 

i (sec.) 




./o 






.1 






•Q) 

o 






-- 






“• !0 













F I ' . i /m \ • T D A < . ' 

C • I ..) ! / \ 1 . . 



- 63 - 










'Kj 

o 

.oi-b. 






Q> 

o 






-.oz 






^.O/l 













o S /o 





T . . 



- 64 - 




/ 0 




iO 


s --- 




^ ^ ^ ^ 






o 


\ 

(rad.) 




-S' 


-to 




- / o 









o 5” fo 




,/o 




. to 


. OS' 




oS 






o 


tJ 

-.05 ^ 


V 

JO^ 




-. or 


‘Qi 

-./o 




to 









f I ' • I ^ ^ > I 

ii' V ,i|' I. . . y 





- 65 - 




. 04 - 




.04. 


.02 vL 




Cl 


r 

6 > 




0 


£>2 


^ * , 


■ — ~Z^02. 






-04 




► 





O 



5 



t C^ec.) 



ro 



/5 





h r. ♦ i ^ r 

V . ,ri 



r • 

f 



I 



- i 



66 - 




O 



s 



i ( sec ) 



S' 





f « I ” T * * . t 

Yl- rsf A .'V \ 



f . 1 . 



-67- 





O 5" /o 'S' 





!- I . 

1 ‘ . 1 * .»!•- 



1 





.2 




. i 




JS*4A ^ C„^* -io 


. / 






— ** 






~.t 


^ S! 

-.j 




-2 









.to 




• to 


? « 

• 05^ 'iL>'i 




.oS 


•<b 

0 ^ 




o 


-•5' 




-OS 


-./o 




~-/o 













69 - 




S.L- 



A/ot I 



iL. 






. o 4 




•4 


.02 




. 02 


05 

o 




• 


02 




I, 

- . 3 a 

A -.02 






- 0 .; 












- 70 - 








.1 




^— . ^ ^ c • • 


. 1 




,A ^ ^ 




W\yW^ 


tJWF^^ 

f \ r X. y" \ 


V.\ / V X 

^\y / 


M “ vr-" 


/ v/ 


-.z 




/ 





r : 

to 




’to 


lb’ \ 






\ 




OS' 


°_z.. V . 




o 


-.05’ 




-.05 


- lO 




“./O 









' I ^ . ' ' ' ' r r j . 

'' ' ~ • ■■ A ■ • , ? T . . 



- 71 - 








04 


.o^l 




.02 


C5 






o 






-.02 




^ * • 30 

oz 


V" S‘ 

-.04 • 




"" o4 


\ 







O S • o 

i (^ec ) 




- 72 - 




to 








- /"\ . 


^ ^ / \ A." 


r 

_ 


N 

\ 

/ 

S 


\y 


0 

-.5- 






-y.o 













O S io ts 










‘db 






0 






*0 y 

-.of "I. 




■S— - - _ X V«. . /' 


~.lo 




/ 










TC 



' L • 



- 73 - 



4^0 

— •— — ■ 




- 


'%jI 0 
s V 






o 






^Zo 






- 4.0 
















I . 



C 






- 74 - 




£/,c sco/e 

ck ® 












- 75 - 




/A# ^co/» 







- 76 - 



appendix I 



Development of the Equations of Motion 

The equations of motion of the airframe are v/ritten 
by equating the forces and moments acting on the air- 
frame to the plane’s reactions in accordance vriLth 
Ne^irton’s lawso The reference axes are considered 
fixed in space. 

The aircraft is considered to be a rigid body; 
aeroelastic effects are not considered. The motion 
is referred to right-handed Cartesian axes. The axes 
used and the sign conventions observed are shown in 
Fig. ^ and Table I. 

Newton’s second law states that the rate of change 
of momentum of a body is proportional to the force 
applied to the body and the rate of change of the 
moment of momentum is proportional to the torque applied 
to the body. The mass of the airplane is considered 
constant for the duration of the dynamic maneuver 
considered. 






1 (^/vv>v) * /WiV" 




u) I L ‘ 

dt 

I M ^ 1- 
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Investigating the dynamics of an infinitesimal 
element of mass drn of the airplane referred to the 
Cartesian coordinate system yields equation set (3)» 
(3) + Pcirr, - iX Rdm - ijX Q dm 

dhy » ( X^) Q dm - X Lj P dm “ 12. dm 

dinj ‘ R dm " 2X Pdm ~ Z Q dm 

For a finite mass these become 
(4-) lij, * P Jcij* + 2’-)dm - ojxqdrn ' R|xzdm 
ky « qJ(Z'*- +X*’) dm - Rjq-?dm — Pj^xdm 

hj - R |(x^ tq9 dm ^ P |zx dm ~ Q jiq dm 

The integral J ( q*--t-£*-) dm is recognized as the 
moment of inertia of the whole airplane about the 
x-axis. Equation set { 4 ) can be written 





PIx - 


Q Ixy 




u * 


Oly - 




- Pixy 


K 




p r„ 


Q Iqt 



The equations of motion relative to inertial 
axes becomes: 

^ 5 i ^ dt 

^ f y = ftVl ^ 

■' dt 

I.fi . rm dW 
dt 

C^) ^L“ ■ P Px ^ ~ ^ ~ x;i ~ ^ Py» ~ P 

2 M ' dhy a, O t Q ly - R 1(^4 ” R.I 

2N1 * dhfc 4-R.ij~P Tyf" "* Q " ^X£tj 

dt 

Equations (6) are written with respect to fixed 
inertial axes. It is expedient to use Eulerian axes 
with their origin at the c.g. which are fixed in the 
airplane and move vrith ito At any instant the airplane 



in a maneuver has a motion referred to the Eulerian 
axes and the axes have a motion referred to fixed 
space. The earth and its atmosphere are regarded as 
fixed space. The absolute acceleration of the aircraft 
referred to the Eulerian axes can be written 



(7) -h -Sx Vx 

Ctt 

where is the total velocity. 

Similarly the absolute rate of change of the 
moment of momentum is written 



( ^ ) d b cu •><. b 

ctt. dt. 

The first term arises from the change relative to 
the Eulerian axes fixed from instant to instant in the 
airplane and the second considers the change of the 
Eulerian axes relative to inertial axes fixed in space. 

Moments and products of inertia measured relative 
to Eulerian axes are independent of time since the mass 
is assumed constant for the duration of a maneuver. 

Many of the terms in equation set (6) go to zero for 
this reason. 

The term of equation (7) can be written; 



(3a) 



to )< ’=■ i j fc. 

P Q R 

U V ^ 

Where t , j ^ k are unit vectors directed along the 

Eulerian X, 2. axes. 

Expanding the cross product and combining with the 
components of the first term of (7) that equation can 



be written in component form as; 

« u + Q w - RV 
v' + R o - Pv^/ 
^ + PV - ao 



( 9 ) 
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By parallel reasoning (6 ) can be written: 



With the Eulerian axes as chosen and shown in 
Fig, B, the xx -plane is a plane of symmetry. There 
is both a positive and a negative value of y for 
each value of i ; consequently and 



The equations of motion of an aircraft referred to 
Eulerian axes can then be written: 



p Fy - V + RIF- PW] 

Z Ft ' +P\/-Q(jJ 

IL- p - R Ix-t + QR (li -Iv) - PQlxt 

ZM - Qly ^ PR(Ix- U) - R'’ 

ZnJ ' RI-J.-P Ixh + P<3 (lt( - fv) lyt 



The left sides of the equations in set (11) are 
the summations of the external forces and moments 
applied to the airplane in flight. These external 
forces are gravity forces, aerodynamic forces, and 
thrust forces. 

The gravity force may be considered to act at the 




- R.' Ivt P‘‘ lyfr -P PQ 
-PR -f 




( 11 ) 




center of gravity making no contribution to the 
summation of moments. 
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The resolution of the gravity force into components 
referred to disturbed Eulerian axes is shovm in Figo 1® 
The force relations may be written 
(12) ^vi[u + Qw - RV] - X 3 

I ^ [V RU - PW] - Yj 
iFi' * +-PV- Qu] - 

by transferring the gravity terras to the right-hand 
side of the equation® The gravity terms Xj, V 3 , and 
are expanded in Figo 1® 

The left side of all the force and moment equations 
now include aerodynamic and thrust forces as well as 
moments due to control surface deflection® 

The disturbed motion of an airplane at any instant 
can be considered the result of disturbing the airplane 
from some steady flight condition® This motion may 
be referred to as disturbed flight and the Eulerian 
axes as the disturbed axes® The forces acting on the 
airplane rUnder these conditions may be regarded as 
the sum of the force required to produce steady flight 
plus an increment required for the disturbance® Accord- 
ingly, each of the total instantaneous velocity com- 
ponents can be written as the sum of a constant velocity 
component during steady flight and a change caused by 
the disturbance: 

{13 ) U = Uc f LA 

V ‘ Vc -<• AT 



It 



«l mi mm* Cii^ 

mm »« «MO. Pf?iy9%i 

§■!««« Jl IViMtfatlllii MfiW^ W 

♦■ j:t“ *•• I 

.'■I’T- ■;%•:.•>- .rj 

,!■ 1 .C- 1 *' - fl| 

■n • I •»««• *11 m»'r*A.mn 4 »l 

V 4tTMi V^^**** ♦W»|y— iiiJ^tJCl - 414^ 

44 4if^ #4 tMliPB «i» 

^mm Im tm* •Q* 

H40IMK4I4I «««lNAt4 I Wgi |wi — 

YMjpmI |«u» 3l» Un»t«>t4| 

40 4^'43j«4i£ II T3u»p*T fM* I iv«liiilfea4 •« aaj 
•>#« 4 4T *>4fUM 'Mhi» «m«1 

«4r^l^ 4«» T'rt^lr 44 of /Li^m<W^ «4 

>o44f#ii fC»* «»«^4fl onC 4i «JH 4B 4v«4 

m mktt»999 t4 ipt priiiliAMHif Hm4I iIAMI 
tAilfl moAM 9»nsa0i^ mi9l 4i# 1^ ^ mfi 

^*4^^ l^tiiir liPlii im ^ ^mmrnk m Mt 

• 4 B« tstm »99 i»i<u fni m «\VSP<- 

«||«M1«« tiu%M*«4 • lb mmm mu mm »ffM4ai| 

<» v>M« • im louiB 44 r># Pm^crmn 

ryi 

- «i * • 







P* Po ^ 1= 

O' 00+^^ 



R. ' i- 



Substituting in (12) considering that s -IR . q 



etc,, and using the expansion of the gravity terras 
from Figo 1 after setting terms like sin 4 ^ »ir> o , 
5in4;-v^ and cos6-i , the set of equations belov; 
can be written. The assumptions that the sine of an 
angle can be replaced by the angle, the cosine is 
equal to one and the product of sines can be replaced 
by zero produces results viiich, although approximate, 
are in good agreement v/ith physical results according 
to Ref, 7o These assumptions, of course, apply to the 
disturbed angles, 

(14) ^ /W1 w. -p Qo W<, Wo Qo^ ” “Vo ft os<^^ TfL 



IL' "p Lx ' /l Lj. CqoKo f 

" (PoQe i- ^<3,-p+ t<l, ) Txfc 

^ P, n. +-R,> - It) Ixt 

+ ZPo-R "P**) 

* /L •) t p Lxt ^ (. Po Oo +Q,')3 't"P^)^1;)"Ix) 

The aerodynamic forces and moments acting on an 



dt dt 





airplane, the left side of the above equations, are 
functions of the flight conditions and of the deflection 
of the control surfaces. 



4 I 1*^ pTi*H||WVhi# f<CI| r*l 9i(i^4aMj»«4r«iS 
«Pn#fl»fiAf — ktwm^rn ^ U» iMt 

. .tM. . fXt^ ) llITVO 

M>.l ••<*», tiU (Ml* k-*-« 

<• IV -mf !■• 4W)«p4aaM« «t Wt 

• I p$tmf4 «uf ^ ^ « *«iMf m iP4 

f 

«tiMv«vf«p 9yamMl» »«jr« # ^ 

<4tn»Mnrwp ij 

^ <4 fiilt W1HU«« #- lll^ 4% 

fA I 4 «m)* 

•• f*t«tf«^3Llk**i» tip* • •! "44yCP 

•I-*- }• A* *y^ 

• 441 ,#«||f lMI# •• - W I II* «^«l *iJ iTfl 1 1 



1 # 
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Each of the forces {/,s>na2. ) and the moments 
{ ) can be expressed for small oscillations 

by expanding by MacLaurin’s theorem neglecting second 
order terms to obtain: 





y y + u. 









where > IK etc* are partial derivatives of 

with respect to etc. Exactly analogous 

expansions are made for V, £, L, and . These 
forces and moments can be expressed: 



X 



(15) K - 
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dir 


^ AV 
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c)U 




t>tU" 


ur 




-t- Ab V 


4- bL 




b/i. 


3\r 


Op 





IX 

^ ur 



ur 



4- ^ s„ 






i>S. 



fU 



M ■' Mo +- 3_M u. «i + aM uj- j_ M 1 • 

N = No + 



M S ‘ 



"a-r 'ait ti? 

Each of the terms in (15) has a physical sig- 
nificance. /„,Yl^and ^ L^^M.^andN. are the forces 
and moments acting along and about the y ^ h , a"d £ 
axes while the plane is in steady flight. The terms 
similar to ^ express the change in the given 

iu. 

force or moment caused by the given disturbance quan- 



tity, The coefficients of the form are the 
basic dimensional stability derivatives. A complete 
discussion of stability derivatives may be found in 
Refs. (1), (2), (3), (7), (B), and (9). 
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Thrust effects are part of the left sides of 
equation set (14) o The equations for steady flight 
condition become: 



(16) y. = t; co^ I 

^ - X 

M, ^ XSj 

Since the Eulerian axes remain fixed vath reference 
to the airplane during a disturbance, the thrust 
components relative to disturbed axes become: 



(17) 



X - T. cos; 

M - T.ij 



where 

AT 

(la) 

and 



T, the thrust during the disturbance, is X> 
can be considered 
aT^ ITlc 












(19) 



X’" X C.OS? (.“=’5)iXu ^ 






^ iVr- 

M = I o ^ 2 j iJ" 2:j m S 

The equations for steady flight are found by 
substituting the expressions obtained above for the 
thrust and the aerodynamic forces (and moments) into 
the left sides of equation set (14) and setting the 
disturbed quantities equal to zero: 

(20) Xg ~ vv-sin ©£ ■+" Xo Cos; “ o 

Yo+ o i- O ^ O 

+ Wcc^Go " Tsct,; - o 



A X « 
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Lo + o o •- o 

Mo 1- o + - o 

k/, r o + o - o 



If now the complete equations for the disturbed 
airplane are written with both steady flight and 

1 

disturbed terras they xvill be of the forms 
Aerodynamic forces + thrust forces = inertia force ■= 
gravity force 

Equation set (20) can then be subtracted from the 
complete equations of motion for the disturbed air- 
plane yielding: 

[] (1 t < 3 o Wo vVo + Qo w + ~ Vj, n- “ v~n. 

~(5 CjovOo fo) 4^ i* C9 cun 6e cm 4c) 

'■ PoVV» -Pour - WoT? " wrl? ~ (.9 'I' ” 3 O=&0 5 w> 4 o " (9^®® 4J 

g - - ^ ’ H - f ^ 

= ly^ [^ur + Pa V. +- Po T +'/o1> - GJol 4 - QoU, “ITo*^ 

-•-(.jscM ea)e- -t C^CoA 4>o)4> w] 



l>Tr “ 



bn. 






^8 



bS„ 



~ ~p' m lyt 4- (.QoP-o t-Qo n_ + Po'l^ )(it ^ i- 'l) 
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Dividing tiie force equations by and the moment 

equations by the appropriate moment of inertia yields 

terms of the form: 

I 



^ M. 



and — 

Ix 



Replacing — by Vw. and J_ bL n/ by L, simplifies 

^ ly 

the notation. These quantities are called dimensional 
stability derivatives, A list of these stability 
derivatives in terms of the dimensionless forms appears 
as Tables III and IV. The flight path angle Vo replaces 
Oo - The equations in final form can be written: 
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APPENDIX II 

Adaptation of the Equations of Motion 
to the Analogue Computer 



The dynamic equations of motion of the airplane 
were developed in Appendix lo These are summarized 



below: 

(c-0 XuU. r 

-^agur +-ur<^ - tIoVo - R,oV* - v/q fv - 

Ya./^ -^VvV 4-Vlp-p 
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These equations are general and allow initial 
velocities, angular or linear, with respect to any 
axiSo Additional control terms, such as flap 
deflection may be included by adding the product of 
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the desired deflection and the appropriate stability 
derivativeo Any attitude of the airplane is also per- 
mittedo 

In order to fit the two analogue computers avail- 
able certain assumptions were necessary. The plane was 
assumed to be in horizontal level flight with no initial 
velocities except IX , forward velocity along the 
X -axis. Stability axes are chosen with the undis- 
turbed X-~axis alligned with the relative wind and the 
y -axis a principal axis. The product of inertia will 
be different from zero. The terms which describe power 
effects are permitted to go to zero by assuming gliding 
flight since these power terms are not essential to the 
dynamic effects which are studied. They merely add 
complexity. 

Solving the resultant equations in terms of the 
highest derivative places them in the form required for 
analogue computer solution; 

(rr-l) u • +- '/■ur •*- ''rn. ~ ~ 6 
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The assumptions made simplify the equations 
sufficiently to permit use of the analogue computerso 
The inertial coupling terms remain® These products 
are handled by the servo-multipliers and permit the 
airframe studied to move \*dth six degrees of freedom® 

The dimensional stability derivatives are obtained 
from Table V, substituting these values in the equations 
above yields the analogue computer equations for the 
basic airframe: 

Longitudinal Equations 

{jr- ^ « — Ot 51 w . O -LGH -w -‘t 'VJl. ~ UTijj-. " 3t.7-© 

(jT- 14 ) Ur *■ > 3 o u + II, ~ "Z-lo ur + (ol.i S^_ 

+ hlo -t- u- ~ ■p'V" 

(IL-i 5) -.000 «u.6<e<j^ -. 0 ? 4 wr‘flfe.a 5j •".Ooi84 ilr +-.888lpri, -.05 

Lateral Equations 

(jr-ifa) Tr * '•&! n. — , 115V 30t-p •*■1.1 .4- 

~ folO n- Hr +ST..‘2.4> 

\^-m) Ip - 4-, 1785 A-~* os 4.-VT +- (a7.,S 

^ ,o755'n_ <Ss c^n. 

^ -h.oiisir - .»«^oa,d.o sos-p “ s ,4-t. 

H- ■“ - .03X3 

These equations of motion are entered directly 
in the analogue computers® The dimensional form of 
the stability derivatives permits the use of real time, 
that is, one second of machine time is equal to one 
second of real time® The coefficients, which are the 
dimensional stability derivatives, are entered as 
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potentiometer settings either directly or, if they 
are greater than one, as a fraction which is sub- 
sequently amplified to full value. The potentiometer 
settings along with the parameters they represent 
are presented in Tables VI and VII, 

The stability derivatives which represent the 
change in a variable due to a control movement are 
based on a one radian control displacement. The 
product of tdie control movement and the corresponding 
stability derivative represents the forcing function 
vdiich disturbs the airframe from its level flight 
equilibrium. Since the equations are entered directly 
with no machine scaling there is a one to one corres- 
pondence between machine units of 100 volts and unit 
movement, or in the case of angular control movement, 
unit radian throw. For a lOO control movement a ratio 
is set up which establishes the potentiometer setting 
for a 10® throw, A lesser or greater movement would, 
of course, be represented by a proportionately greater 
or lesser potentiometer setting. An example is calcul- 
ated below to obtain the potentiometer setting represent- 
ing for a 10® control movement, 

7Cd 6 ; ri4e 

QY /V OF I 
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That is to say that when one machine unit, or 
100 volts, is placed across the potentiometer set 



at the analogue of the change to the pitching 

moment caused by an elevator throw of one radian will 
enter the equations of motion as a forcing function,, 
Ten degrees would be this value multiplied by the 
fraction -12L or .134® The other potentiometer 



settings representing the change caused by control 
movements are similarly computed. For the basic air- 
plane these values appear in Tables VI and VII. 

The computer diagram for the simultaneous solution 
of these six non-linear equations of motion is presented 
as Fig. 9® The inertial coupling terms, vdiich are the 
product or moment of inertia factors and velocity terms 
such as ) , are handled by the servo-multi- 



plier units as indicated in Fig. 9® 

A check was made on the analogue computer solutions 
obtained for the basic airplane (Case I) at sea level 
using equations {III-2C) of Ref. 1 for the longitudinal 
motion and equations (III-97) for the lateral motion, 
nr-'z.e , Ret. L 
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